Effects of spin-elastic interactions in frustrated Heisenberg antiferromagnets 
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The Heisenberg antiferromagnet on a compressible triangular lattice in the spin-wave approx- 
imation is considered. It is shown that the interaction between quantum fluctuations and elastic 
degrees of freedom stabilize the low symmetric L-phase with a collinear Neel magnetic ordering. 
Multi-stability in the dependence of the on-site magnetization is on an uniaxial pressure is found. 



There is a long standing interest in properties of frus- 
trated quantum antiferromagnets because they display a 
large variety of behavior: e.g. conventional spin order- 
ing, spin liquid ground states, or chirality transitions [Q. 
The antiferromagnetic Heisenberg model on a triangu- 
lar lattice is the simplest two-dimensional frustrated sys- 
tem. Quasi-two-dimensional triangular lattice antiferro- 
magnets are not rare objects. It was recently shown [|]||] 
that the materials AM{SO±)% belonging to the Yavapai- 
ite family may be considered as realizations of a quasi- 
two-dimensional triangular lattice quantum {e.g. M=Ti 
with S=l/2) and quasiclassical (e.g. M=Fe with S=5/2) 
antiferromagnets. The rhombohedral /3-phase and mon- 
oclinic a-phase of solid oxygen consist of weakly interact- 
ing planes (triangular in the case of (3 — O2 and rectangu- 
lar in the case of a — O2 || ) whose magnetic properties 
are described by the easy-plane Heisenberg antiferromag- 
netic model with spin S = 1 Q . 

Unconventional properties of the k — (BEDT — 
TTF)2 X family can be described by a Hubbard model on 
a triangular lattice with an anisotropic interactions ||. 
An antiferromagnetic Heisenberg model was also used in 
to describe the properties of a triangular monolayer 
of Pb and Sn adatoms on the (111) surface of Ge form. 

From the theoretical point of view the antiferromag- 
netic Heisenberg model on a triangular lattice has at- 
tracted particular attention since Anderson's proposal of 
a resonance- valence-bond state for this model ^ . From 
that time much work has been done to understand the 
nature of its ground state. There is now an almost 
consensus based on variety of methods |9|-[T^| that the 
ground state has a conventional three-sublattice order 
but the situation is close to marginal and quantum fluc- 
tuations arc important |13[| . Quite recently the effect 
of quantum spin fluctuations on the ground-state prop- 
erties of spatially anisotropic Heisenberg antiferromag- 
net in the linear spin-wave approximation was explored 
jig] . The staggered magnetization and the magnon 
dispersion were calculated as functions of the ratio of 
the antiferromagnetic exchange between the second and 
first neighbors, J2/J1. As a physical tool which allows 
to tunc the ratio J2/J1 Merino et al invoke an uni- 
axial stress within a layer whereby the relative distances 
between atoms are changed. 

In this paper we investigate the ground-state proper- 
ties of the compressible frustrated antiferromagnets, i.e. 



those systems where the coupling between magnetic and 
elastic degrees of freedom is taken into account, under the 
action of uniaxial pressure. We treat the system self- 
consistcntly considering the spin subsystem in the linear 
spin-wave approximation and using the continuum ap- 
proach for the elastic degrees of freedom. We show that 
as a result of this method the behavior of the on-site 
magnetization as a function of the pressure differs sig- 
nificantly from the behavior which was obtained in the 
framework of the spatially anisotropic Heisenberg model 
|^5| , ^6| . In particular, the on-site-magnetization does not 
vanish in the whole interval of stability of the collinear 
phase. 

We consider the two-dimensional triangular lattice. 
The Hamiltonian of the spin subsystem, interacting with 
the displacements of the magnetic atoms, is 

H = — JgSfiSfi+a (1) 
n,a 

Here Sn is the spin operator of the atom located in 
the n-th site of the triangular lattice: n — n\C\ + 
^202(^1,^2 = 0, ±1,±2, ..) where c*i = (1,0), c 2 = 
(I' 2 ) are the- basic vectors of the triangular lattice, 

Js = JO- -r)a(un+a - Uft)) = J (l -r)J2i,j a * u ij a j) 
is the exchange integral (J > is the exchange con- 
stant), u.fi are the atom displacements from their equi- 
librium (without magnetic interaction) , and a is the vec- 
tor which connects a site with nearest neighbors, rj = 
— din. J/d\a\. We consider a small spatially smooth de- 
formation. Therefore it was convenient to introduce the 
components of the strain tensor u^: (un+a ~ ^n)i = 

5 Ej (fit + IS) a i = Ej u v a h (hi = x,y). 

The elastic energy of the two-dimensional triangular 
lattice in the continuum approximation coincides with 
the elastic energy of the isotropic plane: $ = ^ K (u xx + 

Uyy) 2 + \l-i ((U xx ~ Uyy) 2 + Aul y )-p x U XX -PyUyy Wl^TO 

K and fi are the elastic modulus of the system and 
Pi (i = x, y) are the component of the uniaxial pressure. 

Let us use the transformation (see e.g. [0) to the lo- 
cal frame of reference S~ = 55 cos an + Si sin an, Si = 
Si cos qn — S~ sin qn, S*H = St in which quantization 
axis for the spins at each site coincide with its classi- 
cal direction (cos qn, 0, singH) which is determined by 



1 



the vector q. By using the Holstein-Primakoff spin- 
representation Si = v^2(4 + M> Si = i^S/2(bt - 
bit), = S - blba with bl, [bo) being the creation 
(destruction) Bose-operator of the spin-excitation on the 
n-th site and applying the Bogolyubov transformation 



bit. 



N -i/2 e -^a (cosh(6» £ ) a £ + sinh(0g) aO for 

these operators (k is the wave vector belonging to the 
first Brillouin zone and N is the number of atoms in the 
crystal), the Hamiltonian ([l]) in the spin- wave approxi- 
mation can be represented in the form H = NE gr + 



where 



E gr = E c i - i ^2 smn 



(2) 



is the ground state energy per atom in the spin- 
wave approximation. The coefficients of the Bo- 
golyubov transformation are given by tanh 2 (0^) = 

(h " 9kf (/fc + SkV 2 wher e 



J k 



S (~2J(q) + J{k-q) + J(k + q) 
2S (j(k) - J(q) 



9l 



with J(k) = h J^a J® C0S ^ ®' an< ^ 

= fk9k 



(3) 



(4) 



is the magnon energy. E c i = S 2 J(q) is the energy of the 
spin subsystem per atom in the classical approximation. 

The on-site magnetization M = (§Z) = S — (fejjbn) can 
be also obtained by using the Bogolyubov transforma- 
tion. As a result we get 



M = S - — V sinh 

N 



(5) 



Our goal now is to minimize the total ground state 
energy of the system 



T = $ 



E = E 



J 



4) 



pa 



(6) 



with respect to the variational parameters which are the 
components of the deformation tensor e\ = r/(u xx — 
u yy )/2, 62 = rju xy and. the vector q. Here the magneto- 
elastic constant k = 4/i/ (j rj 2 ) characterizes the stiff- 
ness of the lattice in terms of the intensity of the spin- 
lattice interaction, p = (p x —p y ) / (j]J) is the dimension- 
less uniaxial pressure. In Eq. (^) we omitted the terms 
of spin-lattice interaction which correspond to the dilata- 
tion (contraction) of the lattice ( (u xx + u yy )/2 ) without 
changing its symmetry. These terms don't change the 
qualitative behavior of the system. The equations for q 
and e, have the form 



oq da 



(7) 



Let us consider first the case zero pressure: p = 0. In 
the classical approximation when S — > oo the energy of 
the system can be represented in the form 



(8) 



^From Eq.(||) we obtain that the set of equations (Q) 
has three types of solutions: 

i) H-phase: qu = 47r/3 (cos <j>, sin </>) , cf) = 0,±7r/3. 
The spin structure corresponding to each qu is a three- 
sublattice antiferromaget. Three different qn represent 
three possible antiferromagnetic domains. The lattice 
structure is an undistorted triangular lattice (ei = e 2 = 
0) with the group symmetry D^. 

ii) L-phase: <Jl = 2n/y3 (sin^, cos(j)) , </> = 0, ±7r/3. 
The spin structure corresponding to each qi, is a two- 
sublattice antiferromagnet. The lattice structure for 
(j) = is determined by the deformation tensors t\ — 
2S 2 /k, 62 — while for two other q^ the deformation 
tensors can be obtained by rotations by ±7r/3. The lat- 
tice is characterized by the group symmetry D^h- 

iii) S-phase: q s = (g s ,27r/\/3) , 

ei = n^ 1 (cosg s + cos(g s /2)) , e 2 = where cos(<7 s /2) = 
( + 32k/S 2 — 5^/8. Two more vectors q and defor- 
mation tensors t\ and e 2 can be obtained from the above 
equation by rotations by ±7r/3. The spin structure is a 
spiral antiferromagnet. 

We consider the stability of obtained solutions for the 
case when 



ei > 0, e 2 =0, < q x < 



2tt 
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In the subspace (^) the stability condition reads 



(9) 



d 2 T d 2 T d 2 T 



d 2 T 



dej d< & \dq x dei 



> 



(10) 



One can obtain from Eqs (J8|) and (10) that 



• the £-phase is stable when k/S 2 < % 

• both H- and L-phases are stable when | < k/S 2 < 
5. For k/S 2 — A the H- and L-phases have the same 
energy while for k/S 2 < 4 the L-phase becomes en- 
ergetically more favorable. 

• the _ff-phase is stable when k/S 2 > 5. 

The spiral S-phase does not correspond to a minimum of 
the ground state energy. In the subspace (^|) the spiral 
phase for | < k/S 2 < 5 corresponds to the saddle point 
which separates two stable phases H and L. 

The wave vectors qn and determine also the mag- 
netic structure of the H- and L-phases in the spin-wave 
approximation. This is not true, however, for the spiral 



2 



phase where quantum fluctuations change the value of 
the spiral wave- vector q s . 
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renormalized magneto-elastic constant: k/S(S+1) 

FIG. 1. The deformation tensor e\ — (u xx — u yy )/2 
for the collinear L-phase (full curve S=l/2, dashed curve S=l, 
dotted-dashed curve S=3/2) as a function of the renormalized 
magneto-elastic constant k/S(S + 1). Multistability for 
the case S=l/2 is shown in the inset. 

Evaluating numerically the integrals in the r.h.s. of 
Eq. (H|) for q = qjj, ei = one can obtain that 
the ground state energy of the symmetric triangu- 
lar phase in the spin-wave approximation is Tr ~ 
J (-§ S{S + 1) + 1.1723 S). We have also found that 
the H-phase is stable for n > kh where 



6.86 , for S = |, 

kh_ _ I 3.85 , for 5=1, 

S 2 ~ 1 3.26 , for S = f , 

2.99 , for S = 2 



(11) 



Considering the low-symmetry L-phase which corre- 
sponds to the wave vector qi = 2n(0, l/v3) we obtain 
from Eqs (||) , (^) that its energy in the spin- wave approx- 
imation is determined by the expression 

^ = -S(S + 1)(1 + 2ex) + + 5 J(ei) 

/(ei) = (2 + £l )4 / dy (1-2^-^ sin 2 y)E(m) (12) 
J 2 + ei 



where E(m) is the elliptic integral of the second kind 

-2 



cos 2 y 



|L9| with the parameter m 

One can obtain from Eq. (|4|) that in the k — > limit the 
magnon dispersion in the L-phase has the form Wi(fc) = 
2 JS (2+ei)y / (5ei - 2)fc2 + 3(2 + ei )k 2 . This means that 

the magnons in the low-symmetry L-phase can exist only 
for ei > 0.4 and only in this interval we may look for ex- 
trema of the function Tl- It is seen from Eqs (|l|) that 
the quantity J\ — J(l + ei/2) (J2 = J(l — ei)) corre- 
sponds to the antiferromagnetic exchange between first 
(second) neighbors in the spin-wave theory of Heisen- 
berg antiferromagnet on an anisotropic triangular lattice 
M,n6[. Therefore the interval 1 > ei > 0.4 corresponds 



to the interval < J2/J1 < 1/2 in the spin- wave theory 
of Heisenberg antiferromagnet on an anisotropic triangu- 
lar lattice. 

In Fig. [I] we plot the deformation tensor t\ which pro- 
vides an extremum of Tl, as a function of the coupling 
parameter k for three spin values S = 1/2,1,3/2. It is 
seen that the equation dT '/ ' de\ = has solution e± only 
for k < kl- The critical value kl which determines the 
interval of the existence of the L-phase depends on the 
spin S. It increases when the spin S increases, e.g. 



S 2 



3.97 , for S = i, 

4.47 , for S = 1, 

4.64 , for S= §, 

4.73 , for S = 2 



(13) 



but even for S — 4 it is still less (kl/S 2 « 4.75) than the 
critical value for this parameter obtained in the classical 
approximation (kl/S 2 — 5). The reason for this is the 
strong contribution to the effective elastic energy of the 
system from quantum fluctuations (the second term in 
the r.h.s. of Eq. (|l^)). The quantum fluctuations change 
both the stiffness of the lattice (the second derivative of 
the free energy ( fl2"| ) with respect to the deformation ei) 
and the constant of the spin-elastic interaction ( the first 
derivative of the free energy). A single- valued monotonic 
dependence 6i(k) is obtained for S > 1. For S — 1/2 the 
dependence becomes multi-valued, i.e. there is an inter- 



val of k where two values and 



S 1 ) 



< e 



(2) 



) of the 



deformation tensor correspond to each value of the cou- 
pling constant k. Only the solution which corresponds to 

(2) 

the deformation e\ corresponds to the minimum of the 
effective elastic energy (|l2|). 

To find the stability region of the collinear L-phase 
we numerically evaluated integrals in the l.h.s. of the 
inequality ([h]) and found that it holds in the interval 
K < kl- Comparing the stability limits for H- and L- 
phases give in Eqs (|ll| ) and ([u|) we find the following in- 
teresting results. For S > 1 we have always kh < «l 
like in the classical approach where the stability regions 
for H- and L-phases overlap. However, for S = 1/2 this 
is not the case since there is an interval kl < k < kh 
where neither the L- nor the H-phase exist. 

The nature of the phase in the interval (kl, kh) or in 
other words for ei < 0.4, cannot be clarified in the frame- 
work of the spin- wave approach. The reason why this ap- 
proach fails is the following. As it was mentioned above 
the vectors q s obtained in the classical approach are not 
solutions of the extrema conditions (Q) in the spin-wave 
approximation. On the other hand, as it is seen from Eqs 
(Q) and (^) the magnon frequency uj-^ is real in the inter- 
val ei < 0.4 only for q obtained in the classical approach. 
For other q values uji is not positive definite. Thus to 
check the stability of the corresponding phase is stable 
one should calculate the magnon dispersion beyond the 
spin- wave approach. 



3 




0.02-1 , 1 , 1 1 1 1 1 1 1 — 

-0.17 -0.07 0.03 0.13 0.23 0.33 

uniaxial pressure p 

FIG. 2. The effect of the uniaxial pressure on the on-site 
magnetization in the L-phase (5 = 1/2, k = 1.315). The 
inset represents the magnetization as a function of the 
ratio J2/J1 = 2 (1 — £i)/(2 + ei) of the effective exchange 
integrals between the second and first neighbors in the 
presence of the uniaxial pressure. 

Let us consider how the L-phase evolves in the pres- 
ence of an uniaxial pressure p. From Eqs (^) and (0) we 
get 

2ir 

g=— (0,1), e 2 = 0, 

p = K e 1 -2S(S+l) + S^pl (14) 

In Fig. U we show results for the on-site magnetiza- 
tion in the L-phase for S — | as a function of the 
uniaxial pressure p. Fig. 2 and in particular its inset 
shows that our results are qualitatively different from 
those of a recent spatially anisotropic triangular lattice 
study |1|[16). In contrast to the (Ji, J 2 ) -model |§(l6| 
where the ratio J%j J\ was a free parameter in our model 
the corresponding parameter in a compressible triangu- 
lar lattice J2/J1 = 2(1 — ei)/(2 + ei) is determined self- 
consistently from the extrema conditions (|7]). It is seen 
that in the interval of the existence of the L-phase we 
have J2/J1 < 1/2 (ei > 0.4 ) and the on-site mag- 
netization is finite. It is interesting that there exists an 
interval where two values of the magnetization M cor- 
respond to each value of the pressure p. One can show, 
however, that only the state with a larger value of M is 
stable. 

We conclude that the Heisenberg antiferromagnet on 
an compressible triangular lattice differs qualitatively 
from the Heisenberg antiferromagnet on an anisotropic 
triangular lattice. Self-consistent treatment of the elastic 
degrees of freedom shows that the spiral magnetic phase 
is unstable in the classical approximation. This approach 
shows also that in the spin- wave approximation the inter- 
action between quantum fluctuations and elastic degrees 



of freedom stabilizes the collinear L-phase. 
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